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_h ' We succeed in reproducing the I = 1 B mesons, f?i(5720), B|(5745), and _B* 2 (5839) that are 

recently reported by DO and CDF, by using our semirelativistic quark potential model, which also 
succeeded in predicting the mass spectra of the narrow D s j as well as broad Dq(0 + ) and D[(l + ) 
particles a couple of years ago. 

Mass of higher excited states {£ — 1, 2) of B and B B mesons, which are not yet observed, is also 
predicted at the first order in p/mj, with internal quark momentum p and the b quark mass mi,. We 
find the corresponding B s j are below BK/B*K threshold and should have narrow decay widths 
' contrary to most other predictions. Also already established states {£ = and £ = 1) of D, D 3 , 

I , B, and B s heavy mesons are simultaneously reproduced in good agreement with experimental data 

■ within one percent of accuracy. To calculate these D/D B and B/B B heavy mesons we use different 
\f} ' values of strong coupling a„. 

o ' 

\Q . PACS numbers: 12.39.Hg, 12.39.Pn, 12.40.Yx, 14.40.Lb, 14.40.Nd 
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I. INTRODUCTION 



Recently discovered were narrow meson states, D s0 (2317) by BaBar and £> s i(2460) by CLEO @, both of which 
were confirmed by Belle |3[. These are identified as j p = + and I + of cs excited (I — 1) bound states, respectively, 
and have caused a revival of study on heavy meson spectroscopy. Subsequently, another set of broad heavy mesons, 
Z?o(2308) and D^(2427) were discovered by Belle collaboration [J]. Those are identified as cq (q = u/d) excited (I = 1) 
bound states and have the same quantum numbers j p = + and f + as D s j, respectively. The decay widths of these 
excited D s j mesons are narrow since the masses are below DK/D*K threshold and hence the dominant decay modes 
violate the isospin invariance, whereas those excited D mesons are broad because of no such restriction as in D s j 
cases. More recent experiments reported by CDF and DO [j| found narrow B and B s states of I = 1, Si (5720), 
£?2(5745), and £?* 2 (5839). These are narrow because these decay through the D-waves. 

There are some discussions [6] that a quark potential model is not appropriate to describe these new states. 
In fact, mass spectra and decay widths of these states could not be reproduced by relativised/relativistic quark 
potential models @, H, In 0, @], for instance, kinetic terms of quarks are taken as positive definite \/p 2 + m 2 
and spin-independent linear-rising confining and short-range Coulomb potentials are taken into account together with 
spin-dependent interaction terms symmetric in Q and q. In even though degeneracy among members of a spin 
doublet in the limit of heavy quark symmetry is taken into account, they obtained masses more than a hundred MeV 
larger than the observed D s j mesons. The nonrelativistic quark potential model applied both to mesons and baryons 
has been so successful in reproducing low lying hadrons (see the review |I0l |) that people were puzzled why it did not 
work for D s j mesons (atom like mesons). 

To understand these heavy meson states, a number of interesting ideas have been proposed so far. One is an effective 
Lagrangian approach. The authors of [TTI.[l2l| have proposed the modified Goldberger-Treiman relation in an effective 
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Lagrangian, by which the mass difference between two spin doublets (0~, 1~) and (0 + , 1 + ) are reproduced. This idea 
is theoretically successful in explaining D s states but in fail for D states though experimentally the relation holds. 
Other approaches proposed to understand these new states are the QCD sum rule [13[ , implications of a two- meson 
molecule |6l [lij . four quark system [HI, EE 113 j lattice calculation [l8| and so on. The Bethe-Salpeter (BS) equation 
can be an alternative way to take relativistic corrections into account. [l|| H(| HH A coupled channel method [22j is 
another candidate to interpret the new heavy mesons though underlying physics is unclear. These approaches are 
not yet completely established even though they may be interesting themselves. Hence, the problem still remains 
unsolved by these approaches and is challenging. See the detailed review (23| on heavy mesons. 

We believe that a quark potential model still powerfully survives if we treat a bound state equation appropriately, 
being able to predict not only mass differences but also absolute values of hadron masses. In fact, long time ago 
before the discovery of _D s o(2317) by BaBar, two of us (T.M. and T.M.) [24[ proposed a new bound state equation 
for atom-like mesons, i.e., heavy mesons composed of a heavy quark Q and a light antiquark q, in which quarks are 
treated as four-spinor Dirac particles, and Hamiltonian, wave functions, and eigenvalues are all expanded in p/rriQ 
so that a light quark is treated as relativistically as possible and nonrelativistic reduction is made for a heavy quark. 
The heavy quark symmetry is taken into account consistently within a potential model. It is remarkable that the 
model could predict the levels of D s o(2317), D s i(2460), Dq(2308), and D' 1 (2427) quite successfully even using rather 
obscure parameter values due to a small number of data (ground states and only a few excited states) at the time of 
publication [24| . 

Recent experiments on B and B s (discovery of states of I = 1, 5i(5720), B|(5745), and B* 2 (5839)) by CDF and 
DO Q serves us a good testing ground whether our formulation, semirelativistic quark potential model, or other 
formulation also works or not. One of the main purposes of this paper is to show that our model can also reproduce 
these B/B s states within one percent of accuracy compared with the experiments. Another aim of this paper is 
to reproduce and predict the whole spectrum of D/TX and B/B s heavy mesons, including higher states, 3 D\ and 
3 Z?2/ 1 -D2, which we failed to have reasonal values in [24j ]. 

In the previous paper [24| , because of the small number of input data in those days, we searched for a local minimum 
which fits with at least the lowest ground states and found that the x 2 analysis led to a negative optimal value of 
b in a scalar potential S(r) = b + r/a 2 . However, in this paper we have obtained a positive b owing to enough 
experimental input data this time using the Minuit. One difference caused by this sign change appears in light quark 
mass, m q = m u = m^, whose effective mass in the lowest order Hamiltonian becomes m q + b. This effective mass 
becomes negative, ~ —38 MeV for b < in 24j and positive, ~ 86 MeV in this paper. Heavy quark masses, m c and 
rrifc, are also affected by this sign change, i.e., we have obtained a few hundred MeV smaller values in this paper than 
those in [24| . Other significant influence of sign change is on the magnitude of the corrections of lower components of 
a heavy quark, which are comparable to upper component contributions in [2~i| (6 < 0), while in this paper the lower 
component contributions are comparable to the second order contributions in p/m,Q as we will see later. These are 
the main difference from [2~ij in the way of data analysis and otherwise we are using the same approach as before. In 
addition to a good result for mass levels of £ = and £ = 1 states of all D, D s , B and B s mesons, a positive b causes 
not only obtaining plausible values for k = 2 (£ = 2) states (corresponding to 3 Di, ^D^^D^) but also obtaining 
relatively small degenerate mass Mq in a spin multiplet so that the corrections become larger in this paper compared 
with [2J|. That is, the sign of b is important to obtain the whole spectra of heavy mesons even though the lower lying 
spectra are correctly predicted in both cases. Another important and different point from the former paper is that we 
treat even the light quark masses (m Ul m^, and m s ) as free parameters in this paper, while they were given as input 
in[24|. 

In this paper, using the refined parameters obtained from recent new experimental data, we elaborate computation 
of the mass spectra of heavy mesons D, D s , B, and B s mesons, including D s0 (2317), D sl (2460), Dq(2308), Z)' 1 (2427), 
£?i(5720), 5^(5745), and B* 2 (5839). We find that in this new analysis all of the masses of not only above mentioned 
new excited (£ = 1) states j p = + , 1 + , and 2 + but also those of other excited (£ — 1) states j p = 1 + and 2 + , and 
ground {I = 0) states j p = 0~ and l - of D, D s , B, and B s mesons can be reproduced well within one percent of 
accuracy. Furthermore, we predict B s q and B s \ are below BK/B*K threshold, which means one can find them as 
narrow states and is the same results as in 25[ although we have about 100 MeV smaller values for + and 1 + states. 
This is contrary to most other quark potential model predictions as in D s j particles. 

The main difference between our treatment of the quark potential model and others is that we have systematically 
expanded interaction terms inp/niQ and that we have taken into account lower components of a heavy quark in a bound 
state. It turns out that systematically expanded interaction terms are quite different from those of nonrelativistic 
and/or relativised models except for [2l[. In [2l| they used the BS equation and projected heavy quark sector on 
positive states, otherwise their final formulation is very close to ours. Another large difference is the magnitude of 
a light quark mass, i.e., other models including (2lj use a constituent quark mass while we got much smaller values, 
e.g., m u ,m d ~ 10 MeV and m s ~ 100 MeV. 

The paper is organized as follows. In the next section, we briefly review our model, where some necessary formulas 
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are presented for readers' understanding and a discussion on a convenient quantum number k is given. In Sec. IIIII 
numerical results are presented in Tables IIIII and IIVI for the mass levels of D and D s mesons including their excited 
states together with the most optimal values of parameters in Table HT1 where we find that calculated masses are in 
good agreement with all existing experimental data. Those of B and B s mesons are also given in the same section, 
Tables fVTl and fVTTl To show smallness of the second order corrections, we calculate D meson mass spectra with the 
second order corrections and give them in Table IVl Section [TV] is devoted to conclusion and discussion, including the 
reason why our semirelativistic quark potential model approach gives good results compared with others. 



II. BRIEF REVIEW OF OUR SEMIRELATIVISTIC MODEL 



Let us start with an effective two-body Hamiltonian for a system of a heavy quark Q with mass tuq and a light 
antiquark q with mass m, [26l I271]. 1 

H = {dq-pq+f3qmq) + (aQ-pQ+f3QmQ)+f3qP Q S(r) 

+ j 1 - \ ["? • "Q + • n)(a Q ■ n)]\ V(r), (I) 

where S and V are a scalar confining potential and a vector one-gluon exchange Coulomb potential with transverse 
interaction and n = f/r is a unit vector along the relative coordinate between Q and q. S and V are given explicitly 
as 

S(r) = ^ + b, V (r) = ~^, (2) 

6 r 

where a and 1/6 are parameters with length dimension and a s is a strong coupling constant. Since a heavy quark Q is 
sufficiently heavier than a light antiquark q which is orbiting Q, it is reasonable to apply the Foldy-Wouthuysen-Tani 
(FWT) transformation [29( to the heavy quark related operators in H of Eq.([T]). One should notice that the FWT 
is introduced so that free kinetic terms of a heavy quark become Pquiq and one can expand the interaction terms 
in terms of p/tuq systematically, 2 which is equivalent to a nonrelativistic reduction of the original model. By doing 
this way, we can get an eigenvalue equation for an atom- like bound state Qq (2~ij 3 

(Hfwt - m Q ) <g> V'fwt = Eip FWT , (3) 

with 



H FWT -m Q = H- 1 +H +H 1 +H 2 + --- , (4) 

where E — E — mQ (E being the bound state mass of Qq) is the binding energy and a notation ® denotes that 
gamma matrices of a light antiquark is multiplied from left with the wave function while those of a heavy quark from 
right. Hi in Eq.Q denotes the i-th order expanded Hamiltonian in p/rriQ and its explicit expressions are presented 
in [24| and also in the Appendix A. Then, Eq.([3]) can be numerically solved in order by order in p/m,Q according to 
the standard perturbation method by expanding the eigenvalue and the wave function as 

E = E- m Q = E% + El + El + ■ ■ ■ , (5) 

V'FWT = + V>1 + V 2 + ■ ■ ■ , (6) 

where a superscript a represents a set of quantum numbers, a total angular momentum j, its z component m and a 
quantum number k of the spinor operator K defined below in Eq. (fTTj) and a subscript i of Ef and ?/>f stands for the 
i-th order in p/rriQ. 

In practice, it is reasonable to first solve the 0-th order non-trivial equation by variation to get the 0-th order 
eigenvalues and their wave functions and then to estimate the corrections perturbatively in order by order in p/mQ 



1 Miyazawa and Tanaka [28| pointed out that a relativistic bound state equation having this effective Hamiltonian can be derived with a 
boundary condition in relative time that is different from that of the Bethe-Salpeter equation. 

2 There is another way of expanding the Hamiltonian. See |3(J where the Bloch method |31|I is used to expand the effective Hamiltonian 
in p/mQ, which is reproduced in Appendix B. 

3 There is a pioneering work of this approach [32ll . where expanded terms of the Hamiltonian in p/ttlq are treated differently from ours. 
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by evaluating the matrix elements for those terms. The 0-th order eigenvalue Eft and the wave function tpg for the 
Qq bound state is obtained by solving the 0-th order eigenvalue equation, 

H <8> Vo = #o = S 9 ■ Pq + Pq (m, + S(r)) + V(r). (7) 

The wave function which has two spinor indexes and is expressed as 4 x 4 matrix form is explicitly described by 

r = (o *$ m (f)) , (8) 
*Uv - -( u i\ y h), 0) 

Jm r \ iv k (r) y.^ J 1 y ' 

where Vj m being of 2 x 2 form are the angular part of the wave functions and Uk(r) and Ufe(r) are the radial parts. The 
total angular momentum of a heavy meson J is the sum of the total angular momentum of the light quark degrees of 
freedom Si and the heavy quark spin ^Sq: 

j = s e + is Q with s t = l + (10) 

where (= iag I2X2) and L are the 4x4 spin and the orbital angular momentum of a light antiquark, respectively. 
Furthermore, k is the quantum number of the following spinor operator K [24| 

K = -0 q (Z q -L+l), K^ m = k^ m . (11) 

Note that in our approach K is introduced for a two-body bound system of Q q, th ough the same form of the operator 
K can be defined in the case of a single Dirac particle in a central potential |33l | . One can easily show [34[ that the 
operator K 2 is equivalent to Sj and it holds k — ± [se + |) or s# — \k\ — i. k is also related to j as 

j = \k\ or \k\-l (fc^O). (12) 

In order to diagonalize the 0-th order Hamiltonian Ho of Eq.0 in the k representation, it is convenient to introduce 
a spinor representation of Vj m which is given by the unitary transformation (Eq. (|14[) below) of a set of the spherical 
harmonics Y™ 1 and the vector spherical harmonics defined by 

" vW+i) 



Y% = -fiYr", yff = / _!__ V^ i ?ff=-irtxfW, (13) 

where n = f/r. These vector spherical harmonics are nothing but a set of eigenfunctions for a spin-1 particle. Yjm 

(A=L, M, E) are eigenfunctions of J 2 and J z , having the eigenvalues j(j + 1) and m. The parity of these functions 
is assigned as (— 1) J+1 , (— l)- 7 , (— iy +1 for A=L, M, and E, respectively, since Y™ has a parity (— 1) J . The unitary 
transformation is given by 

Y, m \ ( y] + 2 \ „ ( B ■ Y^ 




where the matrix U is defined by 



VWTT \ -Vj V7TT 



Now, substituting Eqs.((8|) and Q into Eq.([7]), one can eliminate the angular part of the wave function, j/f m , from 
the eigenvalue equation and obtain the radial equation as follows, 



•v 



S + V -d r + ± W u fc (r) \ _ _ fc / u fc (r) 



+ f -m g -5 + y; \v k (r) J - E ° [v k (r) ' ' (16) 

which depends on the quantum number A: alone. This is just the same form as a one-body Dirac equation in a central 
potential. Since K commutes with Ho and the states have the same energy Eq for different values of j, these 
states having a total angular momentum j = \k\ — 1 and a spin doublet, are degenerate with the same value of k. 
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TABLE I: States classified by various quantum numbers 



f 


0" 


1" 


0+ 


i+ 


1 + 


2+ 


1" 


2" 


k 


-1 


-1 


1 


i 


-2 


-2 


2 


2 


2s + l; 

h 


1 - 

2 

'So 


i - 

2 

3 Sr 


1 + 

2 

3 Po 


i + 

2 

3 Pi, 'Pi 


3 + 
2 

'Pi, 'Pi 


3 + 
2 

3 P 2 


3- 
2 

3 Cl 


3 - 
2 

3 P>2, 'P»2 


*7 










*~ 2 


*- 2 


* 2 





ft • I 1 • 

()• 



/ n,, -li. .S -li \ • ■ fc = -l I 

no p/mq corrections | ( ^ ^ q s ^ q \ 0" 



\ Chiral Symmetric / ^ Heavy Quark Symmetric J (p/m Q corrections) 

FIG. 1: Procedure how the degeneracy is resolved in our model. 



The parity of the eigenfunction ipQ is defined by the upper ( "large" ) component of Eq. @and hence taking into 
account the intrinsic parity of quark and antiquark, the parity of the meson Qq is given by [34j | 

P= (17) 

which is equal to the parity ire of the light degrees of freedom in the heavy quark effective theory (HQET) as shown 
in Table I. It is remarkable that both the total angular momentum j and the parity P of a heavy meson can be 
simultaneously determined by k alone. Or the states can be completely classified in terms of k and j as shown in 
TableD 

Classification by the quantum number k is convenient to discuss the relation between the symmetry of the system 
and the structure of mass levels. Now, let us look at how masses of Qq mesons are generated in our model. First, 
when the light quark mass mq and the scalar potential S(r) are neglected, the 0-th order Hamiltonian H has both a 
chiral symmetry and a heavy quark symmetry and the eigenvalues with the same value of \k\ are degenerate in this 
limit, where all masses of 0~, 1~, 0+, and 1+ states with \k\ = 1 are degenerate, as m(0~) = m(l~) = m(0 + ) = 
m(l + ) = ?TiQ.[25| Then, by turning on the light quark mass and the scalar potential, niq ^ and S(r) ^ 0, the 
chiral symmetry is broken and the degeneracy between two spin doublets (0~, 1~) and (0 + , 1 + ) is resolved. At this 
stage, members in spin multiplets are still degenerate due to the same quantum number k. Finally, by adding p/m,Q 
corrections, the heavy quark symmetry is broken and the hyperfine splitting occurs as shown in Fig. [TJ 



III. NUMERICAL CALCULATION 



To obtain masses of Qq states, we first solve the eigenvalue equation Eq. p^|) by variation to get the lowest eigenvalues 
and the wave functions. By taking account of the asymptotic behavior at both r — ► and r — > oo, the trial functions 
for Uk{r) and Vk{r) are given by 



Uk(r), Vk(r) ~ Wk(r) exp -(m q + b)r - i (J^J 



(18) 



where 



A = \ k 2 



3 



(19) 



and Wk (r) is a finite series of a polynomial in r 



N-l 



(20) 



i=0 
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TABLE II: Most optimal values of parameters. 



Parameters a c s a" s a (GeV" 1 ) b (GeV) m u , d (GeV) m s (GeV) m c (GeV) m b (GeV) 
0.261 0.393 1.939 0.0749 0.0112 0.0929 1.032 4.639 

which takes different coefficients for Uk{r) and Ufc(r). After getting eigenvalues and wave functions for the 0-th order 
solutions, we use those wave functions to calculate p/rriQ corrections in order by order. Corrections are estimated by 
evaluating matrix elements of corresponding terms in the higher order Hamiltonian. Details of the prescription for 
carrying out numerical calculations and getting reliable solutions are described in [24| . 

In the work of [24(, by fitting the smallest eigenvalues of the Hamiltonian with masses of £>(1867) and D*(2008) 
for cq (q — u,d), D s (1969) and D*(2112) for cs, and B(5279) and 5(5325) for bq, a stron g co upling constant a s and 
other potential parameters, a and b, were determined. At that time following the paper [321 ] . we started the search 
for a set of parameters with b < 0. Using those parameters obtained this way, other mass levels were calculated and 
compared with the experimental data for Di s \/Bi s \ mesons. Light quark (u and d) mass was taken to be lOMeV as 
an input, being close to a value of current quark mass. 

Here in this paper since we have now a plenty of the observed data, we follow a different way of analysis. 

1) First, we determine the six parameters, i.e., a strong coupling constant a c s , potential parameters a and 6, quark 
masses m Uy d 4 , m S) and m c by fitting calculations to the observed six D meson masses, i.e., those of two I = 
(CP, 1~) states and four I = 1 (0 + , 1 + , 1 + , 2 + ) states, and similarly six D s meson masses, using the Minuit \ 2 
analysis. The results are given in Table [TTJ As written in the Introduction, we have found that the optimal 
value for b becomes positive contrary to the former paper 24], which affects all other parameter values as well 
as the calculated values, M ca i c , Mo, Pi, rii, and Cj whose meanings are explained below. 

2) Then, using the optimal parameters obtained this way except for a strong coupling, the bottom quark mass 
nib and a strong coupling a b s which is assumed to be different from a c s for D/D 8 mesons are determined by 
fitting the four observed B meson masses and two B s meson masses. The most optimal values of parameters 
are presented in Table iHl at the first order calculation in p/m,Q. 

The masses with the same value of k degenerate in the 0-th order for members of each spin doublet, are labeled as 
Mq. Ci denotes the z-th order correction in p/rriQ expansion (in this paper n = 1,2) and thus, the calculated heavy 
meson mass M ca i c is given by the sum of Mq and corrections up to the n-th order, 

n 

Afcalc = Mo +J2( Pi+ ^ ' ( 21 ) 

i=l 

where pi includes the contributions stemming from only upper components of a heavy quark, rii from lower, and 
Ci = Pi + n-i in each order. Here one should notice that Mq = tuq + Eq. The degeneracy in k or for the states in the 
same spin doublet is resolved by taking account of higher order corrections in p/rriQ expansion. See Fig. [T] 

In our formalism, each state is uniquely classified by two quantum numbers, k and j, and also approximately 
classified by the upper component of the wave function of Eq.® in terms of the conventional notation 2s+1 £j. In 
terms of this notation, the j p = 1 + state is a mixed state of 3 Pi and 1 P\, while the j p — 2~ state is a mixed state of 
3 D2 and 1 r>2- We approximately denote them with double quotations in Tables llHTlVIII The relations between them 
are given by Q 

/ i»3dl»\ \ 1 / . fo 1 \ / |3e>_\ \ 

(22) 




l" 3 ^") ) = _L v/3 V2 
PW) / V5 I -y/2 V3 




(23) 



4 In the numerical analysis, we set m u = raj in order to reduce the number of parameters. 
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TABLE III: D meson mass spectra (units are in MeV). 
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TABLE IV: D s meson 


mass spectra (units are in MeV). 






Lj(J ) 


Mo 


ci/Ma 


Pi /M 


nx/ Mq 


Alcaic 


Mobs 


l So(0-) 


1900 


0.352 xlO -1 


0.270 xlO -1 


0.816 xlO -2 


1967 


1969 


3 5i(l~) 




1.102 xlO -1 


1.098 xlO -1 


4.076 xlO" 4 


2110 


2112 


3 P)(0+) 


2095 


1.101 xlO~ x 


1.027 xlO~ x 


0.740 xlO" 2 


2325 


2317 


" 3 Pi"(l+) 




1.779 xlO" 1 


1.752 xlO -1 


2.620 xl0~ 3 


2467 


2460 


" 1 Pi"(l+) 


2239 


1.274 xlO" 1 


1.270 xlO" 1 


3.860 xlO" 4 


2525 


2535 


3 P 2 (2+) 




1.467 xlO~ x 


1.466 xlO~ x 


1.035 xlO" 4 


2568 


2572 


3 P!(1-) 


2342 


2.032 xlO -1 


2.008 xlO -1 


2.382 xlO -3 


2817 




" 3 P 2 "(2") 




2.196 xlO" 1 


2.195 xlO" 1 


0.989 xlO" 4 


2856 





A. D and D s 

Calculated and optimally fitted values M ca i c of D and D s meson masses with the first order corrections are presented 
in Tables Hill and llVl respectively, together with experimental data M b s - Those mass spectra of D and D s mesons 
are also plotted in Figs. [5] and [3] and Figure [3] also includes DK/D*K threshold lines. As one can see easily, our 
calculated masses M ca i c are in good agreement with each observed value M b s , and the discrepancies are less than 
1%. Especially it is remarkable that newly observed levels 3 Po(0 + ) and " 3 Pi(l + )" of both D and D s mesons are 
reproduced well, where the masses of D s j mesons are below DK/D*K thresholds. These levels of D and D s mesons 
cannot be reproduced by any other quark potential models which predict about 100 ~ 200 MeV higher than our 
values even though they also succeed in reproducing other levels, i.e., 1 S'o(0~), 3 Si(l~), " 1 Pi(l + )", and 3 P2(2 + )- 5 
This result gives us a great confidence that our framework may give a good solution to the narrow D s j puzzle. 

The first order corrections for D and D s mesons are moderately large, whose amount is an order of 10%. Since 
the charm quark mass in our fitting is rather small, second order corrections might be necessary to get more reliable 
results. However, we have found that they are only an order of 1% or less and thus have neglected them in this work. 
To demonstrate this smallness, we also present the fit for D meson with the second order corrections in Table fVl as 
a typical example. When one looks at this Table, one notices that the fit with the experiments does not increase 
accuracy compared with Table IIIII of the first order. Hence from here on in our computation we use only the first 
order calculations to obtain any levels including radial excitations which we present in a separate paper [351 ]. 

B. B and B s 

For B and B s mesons, unfortunately there are only a few data available although newly discovered states (Pi (5720), 
P|(5745), and P* 2 (5839)) are reported recently by DO and CDF.[5J Calculated masses of B and B s mesons are given 
in Tables PVTl and IVlH respectively. We have used a new value of a strong coupling a b s given by Table ITU for calculating 
B and B s mesons different from a c s for D meson in Table llll Newly discovered levels are well reproduced by our model, 



5 See Table XIV of Hj. 
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too, as shown in these Tables. We predict the mass of several excited states which are not yet observed. The mass 
spectra of B and B s mesons are also plotted in Figs. Q] and [5] and Figure [5] also includes BK/B*K threshold lines. 
Among these levels, predicted masses of + and 1 + states for B s mesons are below BK/B*K threshold the same as 
the case for D s j mesons. Therefore, their decay modes are kinematically forbidden, and the dominant decay modes 
are the pionic decay which is the same results obtained in [25j although values of B s mesons are slightly different from 
each other: 

P s (0+) -> B s {0-)+tt, (24a) 

B,(l+) -> fl.(l-)+7r. (24b) 

The decay widths of these states are expected to be narrow the same as D s j meson cases, since those decay modes 
of Eq. (f2~4"|) violate the isospin invariance. These higher states might be observed in Tevatron/LHC experiments 
or even in DO and/or CDF of FNAL in the near future by analyzing the above decay modes. This is because 
DO and CDF have recently announced the discovery of Pi(5720) = S(" 1 P 1 (1+)"), 52(5745) = P( 3 P 2 (2+)), and 
P* 2 (5839) = B s ( 3 P 2 (2 + )), which have narrow decay width with the same decay products as in Eg. (|2"4")) since only the 
D-wave decay channel is possible. @ Even though these mass values are reproduced by other models [i| 
but again those models predict larger mass values for + and 1 + states of B s mesons than ours, the same situation 
as in D s j so that these states have broad decay width contrary to our prediction. We hope that our prediction could 
be tested in the forthcoming experiments which distinguish ours from other models. 



IV. CONCLUSION AND DISCUSSION 



Given recent new data and expanding a parameter space with positive b, we have reanalyzed the mass spectra of 
heavy mesons, D, D s , B and B s with our semirelativistic quark potential model, in which an effective Hamiltonian has 
both a heavy quark symmetry and a chiral symmetry in a certain limit of parameters. Hamiltonian, wave function, 
and eigenvalue are all consistently expanded in p/rriQ, and a light antiquark is treated as four-spinor Dirac particle 
while a heavy quark is nonrelativistically reduced using the FWT transformation. Calculated masses of not only 
newly observed excited {£ = 1) states, P s0 (2317), P s i(2460), 1^(2308) and L>i(2427) which are not explained by 
conventional quark potential models, but also already established (£ = and I = 1) states of D and D s mesons are 
simultaneously reproduced in good agreement with experimental data within one percent of accuracy. Now we also 
have plausible values of masses for 3 Di and " 3 D2" which we failed to give in the former paper [24|. This has been 
achieved by finding the most optimal parameter set with positive parameter b. 

These results suggest that our model could be a good solution to the narrow D s j puzzle, which is considered to be 
still challenging and people are still looking for exotic state possibilities. We would say that the quark potential model 
still remains powerful enough to understand the heavy meson spectroscopy. The important thing in a heavy-light 
system is to expand all quantities in the model in p/rriQ and set equations consistently order by order and not to 
neglect the lower components of a heavy quark, which has not been taken into account by other quark potential 
models. We have also predicted the masses of higher excited states of B and B s mesons which are not yet observed, at 
the first order in p/m&. Among these, the newly discovered states by DO and CDF are also reproduced by our model, 
too. There are other models to reproduce these states 0, Hi], [3(| [331 but they cannot predict + and 1 + states less 
than BK/B* K threshold. This is the main difference of the results when comparing our model with others. Other 
models for heavy mesons cannot yield the masses of 3 p)(0 + ) and " 3 Pi(l + )" states of D s and B s mesons much lower 
than ,,1 P 1 (1+) V and 3 P 2 (2+) states like ours. 



TABLE V: Second order D meson mass spectra (units are in MeV). 



2s+ % if) 


M 


a /M 


c 2 /M 


M calc 


M obs 


^0(0") 


1783 


0.678 xlO" 1 


-2.091 xlO -2 


1867 


1867 


3 Si(l-) 




1.245 xl0 _1 


2.507 xl0~ 3 


2009 


2008 


3 Po(0+) 


1935 


1.694 xlO" 1 


1.567 xl0~ 2 


2293 


2308 


" 3 Pi"(l + ) 




2.132 xlO" 1 


1.665 xl0~ 3 


2350 


2427 




2045 


1.412 xlO" 1 


2.499 xl0~ 3 


2432 


2420 


3 P 2 (2+) 




1.616 xlO" 1 


-1.003 xlO -3 


2448 


2460 


3 D 1 (1-) 


2127 


1.863 xlO" 1 


2.109 xl0~ 4 


2803 




" 3 L> 2 "(2-) 




2.021 xlO -1 


-1.409 xlO" 3 


2726 
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TABLE VI: B meson mass spectra (units are in MeV). 



2s+1 Lj(J p ) 


Mo 


Pl/M 




Pi /M 




m/Ma 




M calc 


M obs 


^0(0-) 


5277 


-0.161 xlO" 


2 


-3.795 xlO" 


-3 


2.187 xlO" 


-3 


5270 


5279 


3 5i(l") 




0.981 xlO" 


■2 


9.706 xlO" 


-3 


1.107 xlO" 


-4 


5329 


5325 


3 Po(0+) 


5570 


0.401 xlO" 


■2 


1.937 XlO" 


-3 


2.072 xlO" 


-3 


5592 




" 3 Pi"(l+) 




1.412 xlO" 


■2 


1.400 XlO" 


-2 


1.227 XlO" 


-4 


5649 




" 1 Pi"(l + ) 


5660 


1.069 xlO" 


■2 


1.066 XlO" 


-2 


0.289 xlO" 


-4 


5720 


5720 


3 P 2 (2+) 




1.364 xlO" 


■2 


1.364 xlO" 


-2 


2.120 xlO" 


-7 


5737 


5745 


3 ^(1-) 


5736 


2.203 xlO" 


■1 


2.202 xlO" 


-1 


4.583 XlO" 


-5 


6999 




" 3 D 2 "(2") 




1.430 xlO" 


■1 


1.430 xlO" 


-1 


2.092 xlO" 


-7 


6556 





TABLE VII: B 3 meson mass spectra (units are in MeV). 



2s+1 Lj(J p ) 


Mo 


PL /Mo 


Pi /M 




ni/Mo 




M ca lc 


^0(0-) 


5394 


-0.302 xlO" 2 


-0.485 xlO" 


-2 


0.183 xlO" 


-2 


5378 


3 5i(l") 




0.853 xlO" 2 


0.844 xlO" 


■2 


9.249 xlO" 


-5 


5440 


3 Po(0+) 


5598 


0.350 xlO" 2 


0.173 xlO" 


-2 


0.177 xlO 


-2 


5617 


" 3 Pi"(l+) 




1.498 xlO" 2 


1.444 xlO" 


■2 


5.396 XlO" 


-4 


5682 


" 1 Pi"(l+) 


5775 


0.978 xlO" 2 


0.969 xlO" 


■2 


8.257 xlO" 


-5 


5831 


3 P 2 (2+) 




1.263 xlO" 2 


1.261 xlO" 


-2 


2.014 xlO" 


-5 


5847 


3 ^(1-) 


5875 


2.949 xlO" 2 


2.898 xlO" 


■2 


5.104 XlO" 


-4 


6048 


" 3 D 2 "(2-) 




0.564 xlO" 2 


0.562 xlO" 


-2 


1.980 xlO 


-5 


5908 



Some comments are in order. 

(i) As shown in Tables HTlTlVIfl the first order corrections to the 0-th order mass, Ci/Mq, are roughly an order of 
10% for D and D s mesons, while 1% for B and B s mesons. Therefore, p/rriQ expansion works much better for B and 
B s mesons. This is quite natural, since the bottom quark mass is heavier than the charm quark mass. In addition, 
it is interesting that the first order correction becomes larger for higher excited states. It means that the relativistic 
correction must be more important in higher excited states, since the higher order corrections in p/rriQ expansion 
represent the relativistic effects of a heavy quark. This is consistent with a naive picture that the inner motion of 
light as well as heavy quarks may be larger in excited states, and therefore the effects of lower components of a heavy 
quark may become much more important to be included. 
Jjii) The mass difference between members of a spin doublet satisfies the following relation (given by Eq.(53) of 

m c (M D *-M D ) = m b (M B *-M B ), (25) 
m c {M Dt -M De ) = m b (M B: -M B .), (26) 

in the first order of calculation. This is because each of two states having the same value of the quantum number k has 
an equal mass at the 0-th order and hence the mass is degenerate at this stage as shown in Fig. [T] The splitting occurs 
by including the first order (p/itiq) corrections originated from Hi in Eq.Q, which is proportional to f k (m q )/mQ 
with the same functional form f k (m q ) for two states, and therefore Eqs.(|25l |26|) are exactly satisfied in our model. 
Similar relations hold for higher spin states with the same value of k. 

(iii) As shown in Table HT1 the light quark masses m Ut d and m s are considerably small, m Ut d ~ 10 MeV and m s ~ 90 
MeV, compared with the constituent quark masses which are adopted in conventional potential models Q • It should 
be noted that these values are not input but the outcome of this analysis, though in theprevious analysis [24| . 
fn Ut d = 10 MeV is taken as an input. This result is consistent with those obtained in [38l l39l |40||. 

This situation might be described as follows. A light quark is moving around a heavy quark in a heavy meson 
wearing small gluon clouds with large momenta in our treatment of quarks so that light quarks have so-called current 
quark mass. This holds only in heavy-light systems and in the other systems, e.g., quarkonium and ordinary hadrons, 
we should use quarks wearing more gluon clouds, i.e., constituent quarks as adopted in usual quark potential models. 
There may be another interpretation bridging our model with light mass with a constituent quark model. 

(iv) Repeatedly saying, the main difference between other conventional quark potential models and ours is that 
our model systematically expands interaction terms in p/rriQ, takes into account lower component contributions of 
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FIG. 2: Plot of D meson masses in GeV. Solid bars and white circles represent calculated and observed masses, respectively. 
Specific values are given in Table Hill 
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FIG. 3: Plot of D s meson masses. The legend is as for Fig. [2] The dashed lines show DK and D'iC thresholds, respectively. 
Specific values are given in Table ITVl 



heavy quarks in the intermediate states when calculating higher order terms, and adopt very light antiquark mass. 
With regard to 1/mg expansion, for instance, the model of [7], [8| has kinetic terms, \/p 2 + to 2 , while they adopted 
nonrelativistic interaction terms which are symmetric in Q and q or in toq and rriq. The model of [9( somehow 
manipulated to hold heavy quark symmetry and expanded interactions in 1/ uiq but adopted constituent masses for 
light quarks, m u , m^, and m s . Even so, their numerical results give about a hundred MeV larger values for D s j 
particles. Both models have not taken into account the lower components of a heavy quark which naturally appear 
in our semirelativistic potential model. 

(v) We have been successful in reproducing experimental data for both heavy mesons D/D s and B/B s . However 
we have also left unresolved problems in our semirelativistic model. When one looks at Table [TTJ one notices that 
values of a strong coupling used for D/D s and B/B s are largely different from each other. Especially that for D/D s 
is too small to be a strong coupling. Also we have not used a running strong coupling although most of other quark 
potential models take into account this effect in some way. 
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FIG. 4: Plot of B meson masses. The legend is as for Fig. [2] Specific values are given in Table IVT1 
Calculated mass values for 3 _Di and " 3 D2" are too large to put them in the same figure with others. 
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FIG. 5: Plot of B s meson masses. The legend is as for Fig. [2] The dashed lines show BK and B*K thresholds, respectively. 
Specific values are given in Table IVTT1 



APPENDIX A: SCHRODINGER EQUATION WITH AN EFFECTIVE HAMILTONIAN 



The original Schrodinger equation is actually modified by the FWT transformation in our framework in order to 
make non-relativistic reduction of a heavy quark. Now we recall the modified Schrodinger equation of Eq. ((3|) , 



(i?FWT -m Q )(g> ?/>FWT = Eljj FWT , 



(Al) 



The Hamiltonian and the wave function are transformed by the following FWT and charge conjugate operator on the 
heavy quark sector: 



where given are 



Hfwt = U c U fwt (p'q)HU F w T (p q )U c \ 
ipFWT = U c U FW t(pq)iP, 

Ufwt(p) = exp (w(p)j Q -p^j = cosW(p) +j Q -pshiW(p), 



(A2) 
(A3) 

(A4) 
(A5) 
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and some kinematical variables are denned as 



p = ^, tanW(p) = E —, E=Jp 2 +m Q , (A6) 

p mq + h v ^ 

P = P q = ~Pq, P=P q = -PQ, Q = P ~P (A7) 

Here we are considering a scattering of a heavy quark Q and a light antiquark q. Hence, p i; p^, and q given above 
represent momenta of incoming, outgoing quarks, and gluon which is exchanged between quarks, respectively. Note 
that the argument of the FWT transformation Ufwt operating on a Hamiltonian from left is different from the 
right-operating one, since an outgoing momentum p'q is different from an incoming one pq. However, in our study 
we work in a configuration space in which momenta are nothing but the derivative operators. When we write them 
differently, for instance as pq and p'q, their expressions are reminders of their momentum representation. Therefore, 
although the arguments of Ufwt and Up WT look different, pq and p'q are expressed by the same derivative operator 

— iV. The momentum transfer q operates only on potentials and provides nonvanishing results. 

The charge conjugation operator U c is introduced to make the wave function U c ip a true bi-spinor, i.e., gamma 
matrices of a light antiquark are multiplied from left while those of a heavy quark from right, which is expressed by 
using a notation ®. 

Then, the Hamiltonian in the modified Schrodinger equation given by Eq. ([3]) is expanded in powers of p/tuq: 

H FWT -m Q = H- 1 +H +H 1 +H 2 + --- . (A8) 
Hi stands for the i-th order expanded Hamiltonian whose explicit forms are given by 

H-x = -(l+/3 Q )m Q , (A9a) 



H = a q -p + f3 q m q - f3 q (3 Q S + jl + ^ [a q ■ a Q + (a, • n)(otq ■ n)]\ V, 

Hi = -t^—PqP + —f3 q a Q .fp+±-q)s+ ■ W 

2m Q m Q \ 2 J 2m Q v 

1 



(A9b) 



2mQ 



'V.) [P+^q] + >q \ v.s,. 



[a q + {a q ■ n)n] V, (A9c) 



h, = (p + \q) *s-^qxp. w Q %s - J^fv -^Kf s Q y 

— {(P + q){dQ ■p)+p[a Q - (p + q)]+iq x pj Q } ■ [a q + (a q ■ n)n] V, (A9d) 



8m Q 



More details of matrix elements of each order in 1 /niq and properties of wave functions with several operators are 
given and evaluated in Ref. 24]. 

Eq.(TA9i} means that the lowest order equation gives (1 + /3q)V> = 0, which implies that a heavy quark is regarded 
as a static source of color so that projection operator becomes (1 ± (3q)/2, the one at rest system. That is, the solution 
given by Eq.Q includes only upper or positive components with regard to a heavy quark. In this sense, upper/lower 
component of a heavy quark has the same meaning as positive/negative energy component. 



APPENDIX B: BLOCH METHOD 



This Appendix is based on the paper [30(, by which we explain Bloch's method [3l| for a degenerate system. We 
will see the Foldy-Wouthuysen-Tani [2{| transformation is automatically carried out and will obtain the systematic 
expansion of the heavy- light system in 1/mq. To distinguish this formulation from the one we have used, we use 
different notations in this Appendix. 
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We expand the eigenvalues of the Hamiltonian with respect to the inverse of tuq. First we divide the Hamiltonian 
as TL = /3qtoq + V, and we treat V as a perturbation, where actually TL = H. The eigenvalues of the unperturbed 
Hamiltonian are simply ±toq, and the eigenstates of TL are almost degenerate. Then we apply the perturbation 
method of Bloch [3l|. In general, let us consider the Hamiltonian TL = Hq + V. We denote the projection operator to 
the subspace M with an eigenvalue E of H as P, and the projection operator to the subspace M with eigenvalues {£j} 
of TL as V, where all £j converge to E in the vanishing limit of the interaction V. Further we choose appropriate vectors 
{i/jj} in the subspace M such that their projections on M. become the eigenvectors of TL, that is, TCPipj = SjVipj. If 
we define an operator hi by VP = UPVP , then we get the following eigenvalue equation, 

H eff (PV^) = £ j (PVtP j ), (Bl) 

where the effective Hamiltonian TL e f f is defined by TL e ff — PTLhi. Here we must notice that this effective Hamiltonian 
is not hermitian. The eigenfunction of TL, Vipj, is obtained from the eigenfunction of Tieff given by (f>j = PVtjjj, by 
multiplying hi from left, that is, 

Vi(j)j = UPVipj = Vipj. (B2) 

Note that Pipj = ipj. Again we must notice that hi is not unitary, therefore, even if 4>j is normalized, lA<f>j is not. 
V and hi are formally expanded as follows. 

oc 

v = p - J2 Sk i VSk - ■ ■ ■ Sk ^ VSk ^ ' ( B3 ) 

rn=l {ki} 

where S^s are defined by 

S Q = -P and S k = % (k = l,2,--- ) with Q = 1 - P, a = E-H . (B4) 
and the sum over {ki} extends over non-negative integers satisfying the condition X^I^i 1 ^ = m - 

oo 

u = Y, J2' Sk i VSk - ■■■ Sk - VR ( B5 ) 

m=l {ki} 

where the sum over {ki} extends over non-negative integers satisfying the conditions XX*=i ki = m and 

Si=i ki > p [p = 1, 2,- • -m — 1). Making use of the expansion formula for hi, we obtain the effective Hamiltonian 

including the second order corrections in 1/a ~ I/tuq as follows: 

Tieff = EP + PVP + PV-VP + PV—V—VP - PvHtVPVP +■■■ . (B6) 

a a a a 2 

where P = (1 + 0q)/2 and Q = (1 - (3q)/2. 
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